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The nature of the dark matter that binds galaxies remains an open question. The favored 
candidate has been so far the neutralino. This massive species with evanescent interactions is now 
in difficulty. It would actually collapse in dense clumps and would therefore play havoc with the 
matter it is supposed to shepherd. We focus here on a massive and non-interacting complex scalar 
field as an alternate option to the astronomical missing mass. We investigate the classical solutions 
that describe the Bose condensate of such a field in gravitational interaction with matter. This 
simplistic model accounts quite well for the dark matter inside low-luminosity spirals whereas the 
agreement lessens for the brightest objects where baryons dominate. A scalar mass m ~ 0.4 to 
1.6 x 10 -23 eV is derived when both high and low-luminosity spirals are fitted at the same time. 
Comparison with astronomical observations is made quantitative through a chi-squared analysis. 
We conclude that scalar fields offer a promising direction worth being explored. 



I. INTRODUCTION. 

The observations of the Cosmic Microwave Background anisotropics [1] point towards a flat universe. The deter- 
mination of the relation between the distance of luminosity and the redshift of supernovac SNela [2] strongly favors 
the existence of a cosmological constant which contributes a fraction fl\ ~ 0.7 to the closure density. The pressure- 
to-density ratio to of that fluid is negative with a value of w = — 1 in the case of an exact cosmological constant. 
Alternatively, this component could be in the form of dark energy - the so-called quintessence - whose simplest 
incarnation is a neutral scalar field $ with the Lagrangian density 



C = ^g^d^dv® -V(9) . (1) 
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Should the metric be flat and the field homogeneous, the energy density may be expressed as 

p=n = ^- + v(<f) , (2) 



whereas the pressure obtains from TJy = —g%j P so that 

<j>2 

P = 

2 



V(*) ■ (3) 



If the kinetic term is negligible with respect to the contribution of the potential, a pure cosmological constant - 
u> = — 1 - is recovered. Cosmological scenarios with quintessence in the form of a scalar field have been investigated 
[3] with various potentials and their relevance to structure formation has been discussed. 

On the other hand, matter contributes a fraction f?M ~ 0.3 to the energy balance of the universe. The nature of 
that component is still unresolved insofar as baryons amount only to [4] 

n B h 2 = 0.02 ±0.002 . (4) 

According to the common wisdom, non-baryonic dark matter would be made of neutralinos - a massive species with 
weak interactions that naturally arises in the framework of supersymmetric theories. This approach has given rise to 
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some excitement in the community and many experimental projects have been developed to hunt for these evading 
particles. The general enthusiasm has been recently refreshed when numerical simulations have shown that cold dark 
matter would cluster in very dense and numerous clumps [5] (see however [6] ) . The halo of the Milky Way should 
contain ~ half a thousand satellites with mass in excess of 10 8 M Q while a dozen only of dwarf-spheroidals are seen. 
The clumps would also heat and eventually shred the galactic ridge. More generally, this process would lead to the 
destruction of the disks of spirals. A neutralino cusp would form at the centers of the latter. This is not supported 
by the rotation curves of low surface brightness galaxies that indicate on the contrary the presence of a core with 
constant density. Finally, two-body interactions with halo neutralinos and its associated dynamical friction would 
rapidly disrupt the otherwise observed spinning bar at the center of the Milky Way. 

Neutralinos may be in jeopardy. New candidates are under scrutiny such as particles with self interactions [7]. 
An interesting possibility is based on configurations of the above-mentioned scalar field $ for which the pressure P 
vanishes. An academic example is provided by the exponential potential 



Vm = \ P % cxp{-/3($o-$)} , (5) 



where the parameter (3 is 



/247rG 

13 = \l-ar • (6) 



In a flat and matter-dominated universe, such a field $ would behave just exactly as cold dark matter and would 
contribute a fraction 0$ = p%/pc to the closure density. More generally, the kinetic energy <E> 2 /2 should cancel the 
potential V (<&) in order for the pressure to vanish and for the fluid to mimic the effect of non-relativistic matter. As 
will be discussed in the next section, this is actually the case when $ behaves like an axion and oscillates coherently 
on a time scale much shorter than the typical durations at stake. Alternatively, the field $ could have additional 
degrees of freedom and rotate in the corresponding internal space. The idea that the excess of gravity inside galaxies 
may be due to a classical configuration of some scalar field $ has already drawn some attention. The discussion has 
nevertheless remained at an introductory level. A condensate of massive bosons with repulsive interparticle potential 
has been postulated [8] to suppress the formation of structure on subgalactic scales. The polytropic index of this 
bosonic halo varies from n = 1 at low density up to n — 3 at high density. The stability and annihilation of such a 
system has been mentioned in [9] and a limit on the quartic coupling constant of <j> has been derived. The relevance of 
scalar fields to the structure of galactic haloes and their associated dark matter cannot be seriously addressed without 
comparing the theoretical rotation curves to the observations. An exponential potential - with negative overall sign 
- is shown [10] to lead to flat rotation curves. A massless and non-interacting complex scalar field is thoroughly 
considered in [11]. The self-gravitating structure of the field is calculated. 

We nevertheless feel that these analysis may be improved in several ways. To commence, a negative potential 
does not seem quite appealing. It actually leads to unphysical situations where the scalar field rolls down the hill 
indefinitely and converts the infinite amount of energy stored into V < into kinetic energy. Then, as mentioned in 
[12], the rotation curves are assumed to be flat up to an infinite distance. In both [10] and [11], the Bose condensate 
extends to infinity and the mass of the system diverges linearly in the radius r. The Minkowski metric is no longer 
recovered at large r. On the contrary, space exhibits a small deficit of solid angle. We feel that such a behavior is not 
realistic insofar as the rotation curves of bright spirals are actually found to decrease beyond their optical radius [13]. 
Another strange consequence is that Newton's gravitation does not apply even when the fields are weak. Matching 
the metric with the Robertson- Walker form may also be a problem. Finally, the agreement between the predicted and 
the observed rotation curves is only qualitative and in the case of [10] is based on just a few examples. The goodness 
of that agreement is not assessed from a quantitative point of view. 

This motivated us to reinvestigate more thoroughly the subject. In the next section, we discuss the general conditions 
under which a scalar field may bind galaxies. We show that the same field cannot easily account - at the same time - 
for the dark matter at galactic scales and for the cosmological quintessence. We explain the reasons which have lead 
us to consider the model scrutinized in section III. The corresponding rotations curves are derived in section IV and 
are compared by means of a chi-squared analysis to the universal curves unveiled by [13]. The results are discussed 
in section V and prospects for future investigations are finally suggested. 
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II. CAN A SCALAR FIELD BIND GALAXIES ? 



In the weak field approximation of general relativity - or quasi-Newtonian limit - deviations from the Minkowski 
metric r)^ = diag(l, — 1, — 1, —1) are accounted for by the perturbation hy. v . In the harmonic coordinate gauge where 
it satisfies the condition 

d a hl - ^h a a = , (7) 
this perturbation h^ v is related to the source tensor 

Sfj, v = T^u — — g^ v T a . (8) 

through 

h, u {r) = -4G / "7 v ' (Fr 1 . 



By V ) j3^/ 

r> 

The stress-energy tensor is denoted by T M „ . The gravitational potential - which in the quasi-Newtonian approximation 
is nothing but - is sourced by Sqq. Note that in the case of pressureless matter, one would have 2Soo = 2q 

and the gravitational potential would just be sourced by the energy density. So, the scalar field generates the same 
gravitational potential as an equivalent cold matter component with energy density 2Soo- For this reason, the later 
quantity can be called the effective density. In the simplest case of the neutral scalar field $ of Eq. (1), it reads: 

p eS ee 2S o = 2 (9 $) 2 - 2V . (10) 

Note that unlike the true energy density, p e ff contains no space-derivatives <9;<I>. 

In order to generate the gravitational potential well that is observed inside galaxies, a distribution of dark matter is 
generally introduced in addition to the baryon population. An excess of binding ensues and matter is tied more closely. 
Should the scalar field <& be responsible for the haloes of galaxies, its effective density p c g would play the role of the 
ordinary cold dark matter density, and should in particular be positive. Since the gravitational potential in a galaxy 
is essentially static, we would have to assume - a priori - that the time derivative $ = do® vanishes. The effective 
density p e g would therefore reduce to — 2 V (<&) and the field potential would have to be negative. This is actually the 
solution suggested by [10]. However, even if we are not aware of a principle that strictly forbids a negative potential, 
we would prefer to avoid such an unusual assumption. Nothing would prevent the system from being unstable in that 
case and we will therefore disregard this option. On the other hand, a static field with positive potential V leads to a 
positive Newtonian potential and therefore to repulsion. It would lessen the attraction of ordinary matter and disrupt 
galaxies. This is not completely surprising since the same positive V is proposed to accelerate the expansion of the 
universe as its contents repel each other. This property has led to the hasty conclusion that a scalar field could not 
bind matter inside galaxies. 

Let us assume however that the scalar field $ varies much more rapidly than the system in which it is embedded. 
Our Milky Way rotates in ~ 200 million years. If the field changes on a much shorter timescale, the associated effective 
density would be felt through its time average 

p cS = 2^ 2 ) - 2(V) . (11) 

The field $ may oscillate for instance at the bottom of the potential well. The pulsation of the corresponding vibrations 
is equal to the scalar mass m and may be derived from the curvature of the potential at its minimum in $ = 



m 



,2 _ 



V"(0) . (12) 



Actually the field behaves just like an ordinary harmonic oscillator whenever the pulsation m is much larger than the 
wavevector k of the self-gravitating configuration. This translates into the condition 

2tt 

— << R , (13) 
m 

where R is the typical length on which $ changes appreciably For the haloes of galaxies, R is of order a few kpc. 
Because $ varies in time like exp (—imt), the kinetic and potential energies are related on average by 
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{vm = (t) • (14) 

The effective density may now induce a gravitational attraction insofar as 

p cff = 2 (V) > . (15) 

Allowing the field <I> to vibrate quickly has led to an overall change of sign with respect to the case considered in [10]. 
Notice furthermore that the associated pressure vanishes as a result of (3) and (14) so that the scalar fluid behaves 
just like non-relativistic matter. Such coherent oscillations have already been considered in the literature, in the case 
of the axion in particular - see also the interesting discussion of a quintessence field with a late oscillatory stage in 
[14]. 

Another illustration of a fast evolving field is to make it rotate in some internal space. We may look for configurations 
where the dark energy itself - and not its time-average - is rigorously static. A complex field with a uniformly rotating 
phase features the simplest realization of that idea 

$ = £M -iwt _ (16) 

If the field is non-interacting but has a mass to, the associated effective density obtains from 

P -f = S 00 = 2$t$ _ , (17) 

where the potential U = to 2 $^ $. The pressure of the scalar fluid may be approximated by 

P ~ $t$ -[/($) , (18) 

when the space derivatives of the field arc negligible. This amounts to assume once again that the typical length 
R of the system way exceeds l/u>. Whenever the condition to ~ u> holds, the pressure is vanishingly small and the 
scalar fluid behaves as a non-relativistic component. The associated effective density becomes p c s = uj 2 a 2 (r) with 
no explicit dependence on the time. The complete model will be discussed in the next section where we will consider 
the possibility of a boson-star like system extending over a whole galaxy and playing the role of a dark halo. 

We conclude this section by pointing out the difficulty to have a common explanation for both the local dark 
matter and the cosmological quintessence in terms of a scalar field. An excess of gravitational binding on galactic 
scales requires the condition 

$ 2 >^($) (19) 

to be fulfilled. Conversely, should the overall pressure P be negative to account for a cosmological constant, the 
potential would have to satisfy the inequality 

Y<V($) ■ (20) 

We conclude that the pressure-to-density ratio w must exceed the value of —1/3 in order for both conditions to be 
simultaneously met. Such a range seems to be already excluded by the measurements of supernovae SNela [2]. 

III. THE SELF-GRAVITATING COMPLEX AND MASSIVE SCALAR FIELD. 

Boson stars have been extensively studied in the past - see for instance [15-18] . For clarity, we will briefly summarize 
the main features of self-gravitating bosons, following closely the presentation of [16]. We are interested in the stable 
and bounded configurations of a complex scalar field obeying the action 

S = J \f=g~d 4 x £{$,0„*} = Jv^gd 4 x {g^d^d^ - [/($)} , (21) 

where the potential U is invariant under the global symmetry 
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$ — > $' = e la <P . (22) 

The conservation of the corresponding conserved current is crucial for the stability of the boson star. Real scalar 
fields have no stable bounded configurations. One can show that all spherically symmetric minimum energy solutions 
depend on time only through a rotating phase so that the complex field $ may be expressed as in Eq. (16) - see for 
instance the appendix in [16]. In analogy with the hydrogen atom, such solutions correspond to the energy eigenstates 
(n, I = 0, m = 0). We will see later how the discrete energy levels n of a boson star are associated with different values 
of the rotation parameter u. The parametrization (16) of $ is compatible with a static isotropic metric 

dr 2 = e 2u dt 2 - e 2v {dr 2 + r 2 d6 2 + r 2 sin 2 9 d V 2 } , (23) 

where u and v depend only on the radius r. The Klein-Gordon equation reads 

~2v J // i f„,/ i / i ^ A 1 1 i , ,2— 2u 



e 



r 



a'} + uj 2 e- zu a - U'(a) = . (24) 



The Einstein equations provide two additional independent equations of motion. We can choose for instance 

2v" + v' 2 + — = -8irG e 2v {W + V + U} , (25) 



and 



+ ~(u' +v') = 8 7T G e 2v {W - V - U} , (26) 



.'2 

-r 

r 



where the kinetic and gradient contributions to the field energy density are respectively denoted by 

W = e- 2u ^ and V = °- . (27) 

2 2 y J 

From the assumption that boson stars are non singular configurations, asymptotically flat and of finite energy, severe 
restrictions can be put on the boundary conditions for u, v and a. In order to avoid any singularity at the origin, 
the radial derivatives it', v' and a' must vanish at r = 0. Actually, since an angular deficit at the origin would imply 
an infinite concentration of energy, we infer the complementary condition v(0) = 0. This appears explicitly in the 
Einstein equations when they are written in a slightly different way [18]. Anyhow, in this paper, we will focus on the 
Newtonian regime for which v can be shifted by a constant term without any modification for the (u , a) solution. 

Spacetime is asymptotically Minkowskian if both metric parameters u and v vanish at infinite distance r. More 
precisely, for bounded configurations, one expects that on very large distances the field will appear as a point-like 
mass M and that Schwarzschild's metric will be recovered 



e" 



= GM/2 . (28) 



Let us focus now on the finite energy condition. The total energy can be inferred from the matter and gravitational 
Lagrangian. The latter quantity can be calculated by subtracting to the Einstein-Hilbert action a surface term - as 
is usually done for bounded gravitational objects - so that 

L G = -L / r 2 dre u + v {v 12 + 2u'v'} . (29) 

2 Cr Jo 

The total energy is the sum of the gravitational energy Eq = —Lq and of the matter energy 

/>oo 

E M = / 47rr 2 dr e u + 3v {W + V + U} . (30) 



o 



Using the Einstein equation (25), one can rewrite the total energy in terms of the metric and integrate exactly: 

E M + E g = - G- 1 I"" r 2 dr e u + v {v" + v' 2 + -v' 



r 



G- 1 lim \r 2 v'e u + v X . (31) 
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Inserting the Schwarzschild asymptotic solution (28), one can check that the total mass is the same as the total energy 
so that M = Em + Eq. Following (28), the mass is also the limit of a slighlty different expression 



M= -G- 1 lim {r 2 v'e v ' 2 \ (32) 

By rewriting this limit as an integral over r, and by using the Einstein equation (25), one is able to express the mass 
or total energy in terms of the field energy density 



E M + E G = Anr 2 dr e 5v ' 2 {W + V + U} . (33) 
Jo 

Notice that the gravitational contribution is then contained in the factor exp(5u/2). For bounded objects, the sum 
W + V + U should therefore go to zero faster than r -3 . Another important quantity is the conserved charge associated 
to the U(l) global symmetry, i.e., the number of particles minus antiparticles N 



N 



/ AKr 2 dr^g^ {i - M $t$)} = - / 4irr 2 dr e u + 3w W . (34) 

Jo w Jo 



The simplest realization of this system occurs with a quadratic potential U = m 2 $t <3>. By inspecting the Klein- 
Gordon equation at large radii, one finds that finite energy solutions may exist only if m > uj 1 . Moreover, as soon as 
u and v - respectively u' and v' - are small with respect to unity - respectively 1/r - the field asymptotically behaves 
as 

a oc r- 1 exp{- (m 2 - uo 2 ) 1 ' 2 r} . (35) 

Dimcnsionless equations are obtained by rescaling the field by the Planck mass and the radial coordinate by mT x - 
which is essentially the Compton wavelength of $ 



(j = v2ttG a , r = r m . (36) 

Because of the symmetries of the action, the particle mass m, the rotation velocity oj and the lapse function e u appear 
in the dimcnsionless equatii 
the rescaled lapse function 



in the dimensionless equations only through the particular combination (u>/m) 2 e 2u . It is then convenient to define 



e" 2 « = 4e" 2M ■ (37) 



Asymptotic flatness imposes a relation between (u>/m) and the value of u at infinity 



— = e-^ 00 ) . (38) 
m 



The solutions can be calculated by integrating a simple system of three variables — a, u and v - from zero to infinity. 
For a given ct(0), with the assumption that v(0) = and all first derivatives vanish at the origin, there is only one free 
boundary condition left, namely the value of u(0). Using an overshooting method, one finds a discrete set of values 
u(0) n - with n = 0, oo - such that u converges at infinity with a and v smoothly decreasing towards zero. The 
resulting configurations are the energy eigenstates of the system. The state with minimal energy is characterized by 
the absence of nodes - of spheres where a(r) = - while each n-excited state has got n nodes. 

Since we will assume that bosons play the role of galactic dark matter, we only need to study the Newtonian regime 
in which \u\ and \v\ <C 1. In this limit, the system has got additional symmetries which facilitate the description and 
classification of the exact numerical solutions. The global order of magnitude of u and v depends on the parameter 
£ defined by £ 2 = 1 — ui 2 /m 2 , with ( « 1 corresponding to the Newtonian limit. Indeed, one can show [16] that 



1 In the opposite case m < u>, the field oscillates at large distance like r _1 sin((tj 2 — m 2 ) 1//2 r). It fills the Universe with an 
infinite amount of energy, unless some truncation mechanism is put by hand. This problem arises in particular when m = 
[11], but not for the solutions considered here. 
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u = 0(£ 2 ) and v = 0(£ 2 ) while u+v = 0(£ 4 ). So, at order 0(£ 2 ), u = — v - as usual in the first-order post-Newtonian 
approximation [19] - and the system follows a simple pair of equations 

u" + - u = 2 a 2 , (39) 
r 

a" + \ a' = 2 u a . (40) 
r 

The solutions are therefore left invariant by the following rescaling 

u — ► k u , 

a — > k a , (41) 
f — > k- 1/2 f . 

This means that any configuration is fully described by its number n of nodes and by the field value (t(0) at the 
origin. In other words, Newtonian solutions with the same number of nodes are related among each other through 
the rescaling (41). The invariance of the solution appears more clearly when relations (39) and (40) are expressed in 
terms of the ratios S — a(r)/a(0) and Ui e d = u(r)/a(Q) 



<d+:C = 25 2 , (42) 



S" + - S' = 2 u red S . (43) 

X 



The length parameter is now described by x = m r y/a(0) = f y/a(0). We conclude that once the number n of nodes 
is specified, the value of M re d at the origin is unique. So is the field configuration S(x). The ratio M re d(0) = u(0)/a(0) 
has been computed for the fundamental state, the first excited states and also in the limit where n — > oo. Our results 
are quoted in Table I and are in good agreement with [16]. For each configuration, u/m could have been calculated 



n 


U re d(0) = «(0)/ff(0) 


(l-ui/m)/a(0) 




N/y/a(0) 





-0.91858 


0.97894 


2.4 


1.2 


1 


-1.2099 


0.916 


5.4 


2.7 


2 


-1.3437 


0.892 


8.4 


4.2 


3 


-1.4282 


0.877 


11.4 


5.7 


5 


-1.5370 


0.860 


17.4 


8.7 


10 


-1.6831 


0.839 


32 


16 


oo 


-5. 


4.35 


oo 


oo 



TABLE I. Scaling factors for u, (1 — ui/m), M and JV for the fundamental n — state and a few n-excited states. These 
numbers are applicable only in the Newtonian limit which is reached when all quantities u, v, u, a and (1 — u/m) are small 
with respect to unity. 

from Eq. (38) but in practice u converges very slowly. We obtain much more precision by taking into account the 
asymptotic Schwarzschild expression (28) which implies that 

e u = 1 - fv! + 0(r~ 2 ) . (44) 

Noticing that u' = u' , we find in the Newtonian limit 

1 — — = lim {u — fu} . (45) 
m i — >oo 

We also compute a dimensionless mass parameter 

M= lim{fV} = ^ , (46) 



7 



and a rescalcd particle number 



S = i »*••*- 5^"- (47) 



While (1 — Lojra) scales as cr(0), M and iV scale as yo(6) with factors depending on n that are given in Table I. 

We now focus on rotation curves inside a toy model of galactic halo consisting only of bosonic dark matter. Any 
baryonic contribution from the disk, the bulge or any other component is neglected here. Test particles with circular 
orbits of radius r have rotation speed v 



d_ 

dr 



$ grav = C V . (48) 



So, the rotation curve is given by cy/ru 1 . In Fig. 1, we plot this quantity for the fundamental and the n = 2,4,6 
states as well as for an extremely excited field configuration with n — > oo. We also show the rotation curves associated 
to the usual fitting functions for cold dark matter haloes. We first consider an isothermal distribution with 

p iso cx {a 2 + r 2 }^ 1 , (49) 

where a stands for the core radius. We also feature the case of the cuspy profile towards which n-body numerical 
simulations point [20] 

Pcusp oc r _1 {a + r}~ 2 . (50) 

In this case, a approximately corresponds to the radius where peak velocity is reached. The overall normalization 
and a have been adjusted in order to match the scalar field halo rotation curves. As can be seen on Fig. 1, the curve 
associated to the cold dark matter isothermal distribution (49) becomes flat at large distances as in the large n case. 
It does not exhibit however the peak which all the other curves feature. At the outskirts of the system, the cuspy 
distribution (50) leads to a decrease of the rotation velocity with V 2 cx log(r)/r while, in the n = 0, 2, 4, 6 states, the 
decrease is typically Keplerian. Near the origin, the scalar field configurations yield a core with constant density as in 
the isothermal case. This is in agreement with recent measurements of the rotation curves of low surface brightness 
spirals. Note the small wiggles of the excited configurations. This could provide an explanation for some peculiar 
rotation curves with oscillatory behavior - as proposed by [11] who obtains similar curves in his massless model. 
However we will not consider this possibility in this paper, since we focus in a first stage on universal rotation curves 
of spiral galaxies [13]. 

At this point, we must say a few words about the stability of such self-interacting bosons. Beside the Newtonian 
limit, a first complication arises from the fact that - for a given number N of particles and number n of nodes - there 
are actually several values of the mass M corresponding to static, bounded and spherically symmetric configurations. 
Only the lowest energy state is stable. This phenomenon occurs above a critical particle number N ~ 0.3 (Mp/rn) 2 . In 
the Newtonian regime, we are much below this scale and, for a given set (N, n), there is a unique static configuration. 
A stronger condition is that the gravitating boson system should be stable against fission into free particles so that 
M < Nm. It was shown in [16] that this criteria is fulfilled by all states in the Newtonian limit, even for large values 
of n. In that regime, M tends towards N m when n — > oo at fixed N (in Table I, the precision on M and N is not 
sufficient to see this effect). Stability against fission is only a necessary condition. More generally, one should check 
stability (i) at the classical level under any small perturbation leaving N unchanged and (ii) at the quantum level 
under tunneling from excited states to the fundamental state. The former analysis was performed analytically by [17] 
who concludes on the stability of all excited states - at least in the Newtonian limit. A non-perturbative analysis 
was performed numerically in [21] and yielded the opposite conclusion. However, as clearly stated in [21], these 
simulations were based on the most general perturbations which were in particular allowed to violate the conservation 
of the particle number N. Therefore, the positive result of [17] seems to apply to our situation. As far as the stability 
under tunneling is concerned, we are not aware of any previous result. Notice anyway that, in order to be conservative, 
we will only consider n = fundamental state configurations in what follows. 

IV. COMPARISON WITH OBSERVATIONS. 

In general, reproducing the observed galactic rotation curves amounts in modeling the contribution of many com- 
ponents apart from the halo and the observed luminous disk, like a HI gas, a rotating bar or bulge, etc. In order to 
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n = 



r 



FIG. 1. Rotation curves inside a galactic halo that consists of a pure self-gravitating scalar field. The corresponding boson 
is massive but has no interactions. The fundamental and n — 2, 4, 6 states are featured together with an extremely excited 
field configuration for which n — > oo. Conventional CDM haloes are also presented for comparison with mass density given by 
relations (49) - isothermal distribution - and (50) - cuspy profile. Each curve has an arbitrary core radius and normalization. 
We choose to normalize the five scalar field solutions to a common amplitude at the first maximum. These solutions possess 
?i+l maxima, followed by a decay in r _1//2 - as for any bounded object. The amplitudes of the first inner maxima are 
approximately the same, while the outer ones are bigger. For n — > oo, the last maximum and the r~ 1//2 behavior are rejected 
far outside the figure, at infinity: we only see a quasi-flat region with small oscillations. 

make a strong statement, we will restrict our analysis to the simplest case of spiral galaxies and on distances smaller 
than the optical radius r opt defined as the radius of the sphere encompassing 83% of the luminous matter. Indeed, 
for spiral galaxies and on such distances, the only significant contributions to the total density arise from a stellar 
disk with exponential density profile, plus the unknown dark halo contribution: one can avoid introducing a plethora 
of free parameters describing the other components. On the basis of such considerations, Persic, Salucci and Stel [13] 
- hereafter denoted PSS - performed a detailed statistical study over about 1100 optical and radio rotation curves. 
They rescaled each rotation curve to the same size and amplitude by expressing the radius as r/r opt and the speed 
as V(r)/V(r opt ). The rescaled curves were found to depend only on a single parameter, the luminosity Mi. Galaxies 
were divided in eleven classes depending on their brightness Mj and the authors provided for each group of spirals 
the average rotation curve in the range r < 1.1 r opt - They showed that for non-luminous galaxies, rotation curves are 
increasing near the optical radius while for brighter objects, they tend to become flat or they even slightly decrease. 
This result is remarkable insofar as the dynamical contribution of the luminous disk - known up to a constant bias 
factor (3 - is slightly decreasing at r opt . The immediate conclusion is that faint galaxies are always dominated by their 
halo whereas bright spirals only need a very small contribution from non luminous matter. This amazing one-to-one 
correspondance between the disk and halo core density, depending only on one parameter (the magnitude) , is generally 
called the disk-halo conspiracy. 

The purpose of this article is to investigate whether or not a non-interacting massive scalar field halo may account 
for the universal rotation curves of PSS. To achieve this goal, we must solve once again the Einstein and Klein-Gordon 
equations, adding to the former the contribution from the luminous disk. In order to keep a sherically symmetric 
metric, we will describe the gravitational impact of the disk as if it was spherical. This approximation is reasonable 
provided that the corresponding contribution to the mass budget of the system remains small - which is the case for 



faint spirals. In the opposite case, a much more complicated metric needs to be introduced in order to describe the 
asphcrical distribution of luminous matter. We therefore add to the time time Einstein equation the disk density 
contribution p. Should the disk be alone, the rotation curve would be [13] 

VL,(r) - Vl sk (r opt ) / 1 - 97( : / ^ 2 ", 43 ■ (51) 
{{r/r opt y + 0.78 2 } 

The density p is easily derived from a simple Newtonian calculation 

A^n „ "Kiisk , r, "^diskVdisk -i;2 /_ \ /( r / r opt) fK r,\ 
47rGr P = — 1" 2 = V disk( r opt) 2 ' ( b2 > 



r opt 



where we introduce the dimensionless function 

4.38 it" - 78 5.64 u 122 

{u 2 + 0.78 2 } 143 ~ K + 0.78 2 } 



/(«) = r n'. „-^1.43 " r 9 , n^o ?1 2.43 ■ ( 53 ) 



Note that this profile is only valid in the range 0.04 r op t < r < 2r op t since relation (51) obtains from the fit of a 
more complicated expression involving modified Bessel functions. Below r < 0.04 r opt , we maintain a constant density 
which would otherwise diverge. We have actually checked that the details of the disk mass density near the center do 
not affect our results. 

Introducing the additional mass distribution p of the luminous disk into the pure scalar case discussed in section III 
simply amounts to modify the time-time Einstein relation (25) into 

2v" + v' 2 + ^ = -2e 2v {(l + e~ m ) a 2 + e~ 2v a' 2 + p) , (54) 

where the dimensionless disk density p may be expressed as 

p = — —p . (55) 

The radius f and the scalar field a have been previously defined in relation (36). The other Einstein equation (26) 
as well as the Klein-Gordon relation (24) are not affected. It is worth studying the limit of these equations in the 
Newtonian regime in order to gain intuition on their scaling behavior. As before, the relation u = —v ~ u applies in 
this regime and the system reduces to 

u" + - u = 2 a 2 + p , (56) 
r 

a" + -a' = 2ua . (57) 
r 

Relation (40) is not modified whereas the disk mass density p is introduced in the right hand side term of expression (39) 
to yield Eq. (56). As before, the scale invariance of the solution becomes more obvious when these relations are 
expressed in terms of the quantities S = a(r)/a(0) and w re d = u(r)/a(0). The system (56) and (57) becomes 

2 



<cd + -<1 = 2S + n > ( 58 ) 

S" + - S' = 2 u red S . (59) 

X 

The disk density enters through the dimensionless parameter 

p _4irG p 

K = ^2-— - = — ^2—- . (60) 

<7 Z (0) ui cr (0) 

Should 1Z vanish, the fundamental n = mode solution for S(x) would be uniquely determined. Taking advantage of 
relations (52) and (53) allows to express the disk contribution TZ as 

(61) 
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The general scale-invariant solution S(x) in the presence of a disk depends now on the parameters 



a = x opt = m r opt \/<r(0) and 7 



Vt 




(62) 



which come into play through the disk density 



1Z {a, 7} = "fa 2 f(x/a) 



(63) 



The radial extension x opt of the matter disk relative to the scalar halo is accounted for by the parameter a. Since 
the scalar field generates a rotation curve whose magnitude V 2 scales as <t(0), the quantity 7 measures the dynamical 
impact of the disk relative to that of the halo. These parameters determine therefore the size and the mass of the 
disk with respect to the halo in which it is embedded. The actual scale of the entire system is specified in turns by - 
say - the optical radius r opt and the disk speed Vdisk^opt)- Once the configuration {0,7} is chosen, the behavior of 
the field - given by the scale-invariant solution S a ~,(x) to the system (58) and (59) - is completely determined. The 
full rotation curve which both disk and halo generate may be readily derived. 

With a \ 2 analysis, we obtain likelihood contours in the two-dimensional 2 free parameter space {a, 7} shown on 
Fig. 2. We use only the six less luminous galaxy classes from PSS since they are the most relevant probes of the 
halo contribution and a priori the closest cases to spherical symmetry - brighter galaxies require only a small halo 
contribution, at least inside the optical radius. The minimum \ 2 values are (2.5, 23.6, 20.9, 29.4, 17.8,7.6), and the 
corresponding best-fit rotation curves are shown on Fig. 3. We will not over-interpret the absolute value of the % 2 in 
terms of goodness-of-fit, because we do not know the exact meaning of the data error bars. A careful examination 
shows that the points are not distributed according to their very small errors, at least if the data are to be explained 
by smooth curves - this is visible for instance with the third point in the < Mj >= —20.5 case, which explains why 
the minimal \ 2 is only 29.4. Therefore, there is a hint either that the errors are slightly underestimated, or that the 
data feature small wiggles corresponding e.g. to spiral arms, that should enter into a better modelization of the disk 
density. Anyway, even with the given error bars, the y 2 are already fairly good and our scalar halo model seems to 
fit universal rotation curves at least as well as the toy-model CDM halo used by PSS. We also note that for the last 
two cases, the minimum is strongly degenerate along a line ranging from a = to a = 0.5. So, for the most luminous 
galaxies, the data provide only a lower bound on the halo size, while for the other cases a specific value is preferred. 

For each individual galaxy, one can restore the physical value of the rotation velocity by adjusting the free scaling 
parameter er(0). The mass of the scalar field is given by 



Here we defined the parameter j3 = V| isk (r op t)/V 2 (r opt ), calculated by the code for each value of the relevant param- 
eters {a, 7}. Of course, our model can provide a powerful explanation for galactic rotation curves only if all galaxies 
can be fitted simultaneously with the same value of m, and therefore approximately the same angular velocity lo 
(in principle one could introduce a time-dependant effective mass, with slightly different values at large and small 
redshifts; however, since galaxy rotation curves have never been found to be redshift-dependant, we discard this pos- 
sibility and assume that m is constant in space and time, at least on observable scales). In order to test this idea, it 
would be natural to use data from individual galaxies; but in doing so, one would face back the problems associated 
with large systematic uncertainties, which motivate the PSS approach. In a first step, we plot \ 2 as a function of 
m r op t V(r op t) in Fig. 4 and remark that the faint galaxies have a preferred mass while the brighter ones in the two last 
panels provide only an upper bound. Then, we assume that each synthetic universal rotation curve is associated with 
a unique typical galaxy, with average optical radius and velocity. For each class of magnitude, we perform an average 
on the sub-sample given by PSS - see the tables in their Appendix D - and find respectively (r opt V(r op t) = 710, 1100, 
1900, 2500, 3200 and 4900 kpc km s _1 . Since we do not employ the raw 600 galaxy data, we cannot give an error on 
these numbers, nor can we make a precise prediction for the mass. However, plotting now the \ 2 as a function of to, 
we find that the first four classes are perfectly compatible with a mass m € [4 — 16] x 10 _24 eV while the two others 
push towards the lower-end of this interval with me [4 — 8] x 10 -24 eV. 



2 The PSS data include an error bar also at r = r opt , reflecting the cumulated observational uncertainties at the point chosen 
for rescaling. So, in each \ 2 calculation, we must marginalize over an overall data normalization factor, which means that we 
have three free parameters and 11 — 3 = 8 degrees of freedom. 
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FIG. 2. Two-dimensional likelihood contours in the {a, 7} space, for six galaxy classes ordered by growing magnitude. In 
the four first cases, the crosses show the best-fit models. For the last two cases, the minimum is strongly degenerate along a 
line ranging from a = to a — 0.5. The curves correspond to \ 2 ~ Xmin — 3.2 (resp. 6.2), which would give the 68% (resp. 
95%) allowed region if the experimental errors could be rigorously interpreted as 1-a gaussian errors. 



On Fig. 5, we plot the rotation curves obtained with m = 6 x 10~ 24 eV ~ 5 x 10~ 52 Mp and minimized over 7. 
Although the % 2 values are approximately 1.5 bigger than those of the independent best-fits of Fig. 3, the agreement 
with the data remains quite good. The effect of fixing the mass is to obtain more radical behaviors at r opt : for light 
galaxies, the rotation curves are growing faster near the optical radius while for bright galaxies they are even flatter. 
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FIG. 3. In each panel, the theoretical rotation curve corresponds to the {a, 7} configuration that provides the best fit to 
the data. The dashed lines shows the contribution of the disk, which increases with the galaxy magnitude as for usual CDM 
models. 



For all these models, the total mass is close to M — 5 x 10 10 M Q whereas the ratio of the total radius over the optical 
radius - which respectively encompass 83% of the total and luminous mass - varies between 4 for < Mi >= —18.5 
and 1.5 for < Mj >= -21.2. 
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FIG. 4. The constellation of {a, 7} configurations is featured in the plane (m (r op tVopt) , X 2 )- Faint galaxies have a preferred 
mass while the brighter ones in the two last panels provide only an upper bound. 



V. DISCUSSION AND PROSPECTS. 



In this work, we solved the Einstein and Klein-Gordon equations for a free massive scalar field, in presence of a 
baryonic disk. Using the universal curves of Persic, Salucci and Stel [13], which are based on hundreds of galaxies, we 
showed that a galactic halo consisting in such a Bose condensate could explain fairly well the rotation of low-luminosity 
spiral galaxies. A single value of the mass, of order 6 x 10 _24 cV, is compatible with galaxies of different magnitudes. 
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FIG. 5. In each panel, the average value (r op tVo P t} over the corresponding spiral galaxies has been used to get the scalar 
field mass m = a^/7/(r op t V op t). The best fit configurations with a common value of the mass m = 6 x 10" 24 eV are shown in 
solid lines, together with the observations. Again, the dashed line shows the contribution from the disk only. 



The corresponding Compton wavelength 1/m ~ 10 2 kpc is three orders of magnitudes smaller than the typical size 
of spirals. Indeed, the spatial extension of a self-gravitating field $ is given approximately by l/(y / $(0) m), while the 
square root of the central field value (expressed in Planck units) is comparable with the velocity of orbiting particles 
(in units of c). Since we are dealing with speeds V ~ 100 km s _1 ~ 10~ 3 c, there is really a factor of 10 3 between the 
Compton wavelength and the size of the halo. Note that since cj ~ m, the light scalar field $ rotates in its internal 
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space with a period of ~ 30 yr. 

We conclude that scalar fields could be a nice alternative to CDM halo models. Our positive results concerning 
the rotation curves are strengthened by other astrophysical considerations. First, a scalar field solves naturally the 
dynamical friction issue for bared galaxies. Because it is completely smooth, such an extended field cannot slow down 
the spinning bars observed at the centers of many galaxies, as would a granular CDM medium. Second, following 
Hu, Barkana and Gruzinov [22], an ultra-light scalar field can avoid the excess of small-scale structure predicted by 
CDM simulations near the galactic center. 

Of course, several improvements are needed before concluding that a scalar field is the best galactic dark matter 
candidate on the market. First, it is necessary to extend the comparison to various types of individual galaxy rotation 
curves, with the drawback that more degrees of freedom must be included in realistic modclizations of the baryonic 
components (gas, bulge, ...). This is however the only way to obtain better constraints on to, and to find out whether 
a quartic coupling, not considered in this analysis, improves the model. In some particular cases, it would be worth 
taking into account excited field configurations, which seem to be stable (due to charge conservation) and which 
predict ultra-flat rotation curves far from the core, with small wiggles that may have already been observed (see also 
the claim in [23] concerning possible existence of discrete dynamical classes for spiral galaxy disks). Also, in order 
to get a better view of the rotation curves in the vicinity of the core, especially for bright spiral galaxies, further 
technical ingredients must be passed to the equations, in order to distinguish the spherical symmetry of the halo from 
the quasi two-dimensional distribution of the stars. 

Finally, it would be extremely interesting to plug such a complex light scalar field into a general cosmological 
framework, and study into details the growth of linear perturbations and the formation of non-linear structures. The 
pioneering discussions on such cosmological scenarios [22,24] are very promising and suggest that many interesting 
developments on scalar field dark matter should arise in the next years. 
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